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INPUT -OUTPUT  STABILITY  ANALYSIS  WITH  MAGNETIC 
HYSTERESIS  NON-LINEARITY  -  A  CLASS  OF  MULTIPLIERS 


M.  C.  Safonov*  and  K.  Karimlou* 


Department  of  Electrical  Engineering -Systems 
University  of  Southern  California 
Los  Angeles.  CA  90089-0781 


Abstract 

A  class  of  positive  real  multipliers  is 
obtained  to  establish  frequency  domain  conditions 
for  stability  of  feedback  systems  containing 
ferromagnetic  hysteresis  non-linearity. 

1.  Introduction 

Popov  Criterion  [1]  and  its  extensions,  such 
as  derivation  of  more  general  classes  of 
multipliers  [2],  [3],  consider  non-linear  elements 
that  are  memoryless  and  pass  through  the  origin, 
i.  e,  ,  *(0)  »  0.  An  important  class  of  non- 
linearities,  ferromagnetic  hysteresis,  is 
neither  non-dynamic  nor  passes  through  the 
origin.  Therefore,  to  analyte  the  stability  of 
systems  containing  this  type  of  non-linearity, 
appropriate  modifications  to  Popov's  approach 
should  be  made. 

Published  material  to  address  this  problem 
is  scarce.  Authors  in  [4]  obtained  multipliers  of 
the  form  1+qjw.  q  >  0  for  finite  gain  L^-stability 

of  magnetic  hysteresis  feedback  systems.  This 
paper  generalizes  the  results  of  [4].  The  only 
other  work  known  to  us  that  addresses  the 
problem  of  stability  of  magnetic  hysteresis 
systems  is  by  Lecoq  and  Hopkin  [5],  where  by 
letting  the  derivative  of  their  input  signals  to 
belong  to  exponentially  weighted  L2~  spaces.  they 

obtained  similar  multipliers,  2+qju.  for  bounded 
input -bounded  output  stability  of  h  ^teresis 
systems.  For  stability  of  systems  with  other 
types  of  hysteresis  non-linearity,  such  as 
backlash  and  relay,  see  Yakubovich  [6],  Hsu 
and  Meyer  [7],  Kodama  and  Shirakawa  [8], 

Maeda,  Qceda,  and  Kodama  [9 J,  (10 J. 

In  the  present  paper  we  obtain  a  general 
class  of  positive  real  multipliers  for  the 
stability  analysis  of  feedback  systems  of  the 
type  shows  in  Fig.  (la),  where  N  is  a  ferro¬ 
magnetic  hysteresis  non-linearity  and  H  is  a 
linear  element.  The  analysis  is  done  by  substi¬ 
tution  of  the  model  for  the  hysteresis  by  Chua 
and  Stromsmoe  [11].  Then  the  concepts  of  posi¬ 
tivity  and  passivity  are  utilized  to  deriv  e  frequency 
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domain  conditions  on  the  linear  element  H  for 
finite  gain  L2-stability  of  the  feedback  system. 

II.  Hysteresis  Modeling 

The  model  for  ferromagnetic  hysteresis  of 
Chua  and  Strosmoe  [11]  is  given  by 


*  go[x(t)  -  fo(y(t)] 


where  x(t)  end  y(t)  are  real-valued,  continuous 
input  and  output  signals  of  the  hysteresis  non¬ 
linearity  representing  the  current  i(t)  and  the 
flux  linkage  cp(t )  of  an  inductor  (transformer); 
and  g  and  f  are  strictly  monotonically 
increasing,  differentiable,  onto  functions 
enjoying  the  important  property  of 

g(0)  *  f(0)  *  0 

Equation  (2. 1)  models  the  behavior  of 
ferromagnetic  hysteresis  successfully  and  with 
very  good  accuracy.  It  predicts  the  expansion 
Of  the  area  of  the  hysteresis  loop  with 
increasing  frequency  and  predicts  minor 
hysteresis  loops  such  as  commonly  occur 
when  a  d-c  plus  periodic  input  is  applied. 

o:  Functional  Composition. 
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After  plotting  the  hysteresis  loop  for  * 
convenient  signal,  simple  procedures  ere  given 
in  [11]  to  determine  g  end  f  for  that  loop.  When 
non-linear  functions  g  and  f  are  determined,  they 
can  be  substituted  in  the  model  of  equation  (2.  1) 
to  predict,  with  good  accuracy,  the  hysteresis 
shape  and/or  its  output  for  any  arbitrary  input. 
For  further  detail  and  examples  see  [llj. 


L^-Stable. 


Proof:  See  for  example  [12]. 

Definition  (3.4)  [121;  Let  H:  L?-  L2.  Then 
H  is  positive  iff 


<  Hx.  x  >  >  0  . 


V  x€L„ 


III.  Mathematical  Preliminaries 

Definition  (3.  1)  [121:  Let  H:  L^. 

Then  H  is  passive  iff  there  exists  some  constant 
0€R  such  that 


<  Hx,  x  >T  >  6  V  *€L2e 
V  T6 J 

a 

Definition  (3.  2)  f  121:  Let  H:  L2e- L2#. 

Then  H  is  strictly  passive  iff  there  exists  6  >  0 
and  some  constant  0  €  R  such  that 


<Hx,  x  >T  >  6  [(Xj.ll  ♦  0  Vx€L2e 
VT€I 

a 

Definition  (3.3)  T131:  The  feedback  system 
of  Fig.  (la)  is  said  to  be  finite  gain  L2 -Stable  if 

.)  e,..2.y,.y2€L2  Vu,.u2€L2 

b)  There  exists  constants  and  p2  such  that 

II  «j  II .  II  «2II  •  lly  j  II  •  Ily2!l  <  s  l  ll«J  !!+p2IIu2II 


V  Uj.u2€I.2 


In  the  following  well-known  theorem,  the 
concept  of  passivity  is  used  to  establish  finite 
gain  L2-*tability  of  the  feedback  system  shown 

in  Fig.  (la),  where  N  and  H  are  considered  to 
be  operators  in  the  general  sense. 

Theorem  (3.  1):  Consider  the  feedback 
system  shown  in  Fig.  (la) 


Uj  -  He, 


*2  *  U2  +  Nel 

where  H,N:  l*2c~  L2e*  Assume  that  for  any 
Uj,u2€L2  there  are  solutions  e^,  e2€L2e* 

Suppose  that  there  are  real  constants  v,  6,  and 
«  such  that 


IIh*||t  <  v|l*!!r 

(3.1) 

< Hx,  x  >T  >  a  ||x||* 

(3.  2) 

<x.Nx>t  >  e||Nxl|* 

(3.3) 

Vx€L2€,  V T€ J* 


Under  these  conditions  If 


6  +  t  >  0 


(3.4) 


Note  that  unlike  passivity,  which  is  defined 
on  positivity  is  defined  on  L2. 


Hysteresis  N 


Table  t 
Noutlom 


Symbol 

Flo*  of  real  and  positive  real  numbers. 

*■> 

L. 

Tha  apaca  at  signals  aucb  that  J  |*(t)  jdt 
delete.  -• 

Tha  apaea  o l  bounded  signals . 

The  spec*  of  signals  which  are  square 
intagrable  op  every  bounood  interval 

[O.T]  fMj. 

The  Hilbert  epaet  of  signals  which  ars 
square  integfnble  on  with  inner 

product  <  a,  1  >• 

<i,y> 

J*  y*(t)  *(D  di 

M 

/c  a, a  > 

>x 

aft)  rft 

««|!t 

*'*>T 

*T 

Truncated  «  (151- 

J 

at 

*«  r* ) 

Sot  of  instances  of  time  of  interest. 
Convolution  Algebra  []2]. 

Kaal  part  of  a  complea  quantity. 

Next,  a  positivity  result  for  a  cascade 
combination  of  a  linear  and  a  non-linear 
operator  is  presented. 

Theorem  (3.2)  [13]:  Assume  G:  L2-L2  is 
a  linear  operator,  with  G(jw)  €  Lm ,  which  maps 
x€L2  into  the  element  in  L2,  and  cpt  L2-L2  is 

a  monotone  non -deer easing  non-linear  operator. 
Let  C(juu)  be  given  by  the  Fourier- Sti el jes 
integral 

c(juj)  *  i  -  r  •~imr 


then  the  feedback  system  is  finite  gain 


dV(r) 


(3.4) 


3 


where  V(t)  is  Any  monotone  non- dec  reeling 
function  of  total  variation  less  than  or  equal  to 
unity.  Then  Gcp  is  a  positive  operator  on 

I 

IV.  Main  Results 


Substitution  of  the  model  given  by  equation 
(2.  1}  for  the  hysteresis  non-linearity  N  gives  the 
feedback  system  of  Fig.  (2).  Note  that  although 
the  standard  magnetic  hysteresis  non-linearity, 
as  shown  in  Fig.  (lb),  is  the  plot  of  flux  linkage 
(Q(t))  vs.  current  (i(t))  of  an  inductor  (trans¬ 
former).  but  from  circuit  analysis  point  of  view, 
the  input  and  output  of  the  model  replaced  for  N 
as  shown  in  Fig.  (2)  are  current  through  and 
voltage  across  the  inductor  (transformer).  ^ 

The  class  of  positive  real  multipliers  is 
defined  next. 


Definition  (4.  1):  1ft  is  the  class  of  multipliers 
M(s)  of  the  form: 


m(s)  =  e js+s  )  +  r 

i*i 


•  +  P; 


where:  0  =0  or  1;  *  ,  a.  >  0;  or.  >0,  Vi.  |#1; 
o  o  l  —  i 

p.  >  0,  Vi;  *  <  Min  {p.}. 

1  o  i  criJ  a 

1ft  is  similar  to  the  class  of  multipliers 
studied  in  [2]  and  [3]  in  stability  analysis  of 
standard  Popov  type  feedback  systems  with 
monotone  non-decr easing  non-linearity.  It  is 
shown  that  [2],  M(s)  can  also  be  represented  in 
the  form  ^ 

n  <»+\j 

M(s)  =  i*l - 


k 

n  (s+t|.) 

i*i  * 


where:  l=k  or  k+1;  0<  ‘ 

The  above  representation  implies  that  M(s) 
consists  of  a  combination  of  real  valued  poles 
and  zeroes  in  the  left  half  plane,  where  the  first 
singularity  is  a  zero,  and  the  poles  and  zeroes 
alternate. 


Next,  the  main  stability  result  is  presented. 

Theorem  (4.  1):  Consider  the  feedback 
system  of  Fig.  (2),  where  h(t):  Lj(R+)3,  and 

h(t)€j(.  Assume  that  for  any  Uj,  u2 €L2,  there 
are  solutions  «2» Yj*  ^2e’  **  *or  ,ome 

constant  $  >  0  and  M€7? 


inf  Ra  [M(j«)H(jw)l  »  6  >  0,  V«i>0  (4.1) 

then  V  uj *  *  ttj€  Lj, 


e  voltage  across  an  inductor  (transformer) 
is  proportional  to  the  rate  of  change  of  the  flux 
linkage,  the  constant  of  proportionality  being 
the  number  of  turns. 


Jh(t)  is  impulse  response  of  H(s). 


»>  (i)  •y  «1»«2»y1*y2»y2^-2; 

(ii)  There  exists  constants  and  p2 

such  that 

lUjIUlejIl.  Ikjll.  BtjII.  Ily2ll.  II  y2l! 

<  llujll  +  lluj  If  +p2  llu^ll 

i.  e. ,  finite  gain  L2 -stability, 
b)  Sj,  02^^ y2€Lw  are  continuous,  and  go  to 
zero  as  t 

B 

Proof:  By  taking  the  truncated  inner  - 
products  of  each  branch  of  the  feedback  system 
and  utilizing  Theorem  (3.  2).  positivity  and 
consequently  passivity  of  each  branch  is  proven. 
Then  using  Theorem  (3.  1),  stability  of  the 
system  is  deduced. 

For  detail  of  the  proof,  see  Appendix. 

B 

Upper  and  lower  bounds  of  non-linearities 
g  and  f  can  be  taken  into  consideration  to 
obtain  less  conservative  classes  of  linear 
element  H(s).  As  an  example,  an  upper  bound 
of  g  is  exploited  in  the  following  Corollary: 

Corollary  (4.  1):  Consider  the  feedback 
system  of  Fig.  (2),  where  g€  sector  (o.k), 

<  k,  h(t)€Lj(R+),  h(t)€j(.  Assume  that  for 
any  there  exists  solutions  ej,e2. 

y^.y^L^  If  lor  some  constant  b  >0  and 

Mtlft 

Re  {M(ji»)[H(jui)  +  i]}  =  6>0.  V  1  >0  (4.2) 

then  V  Uj  ,  iij ,  u2*  u2€  conclusions  of  Theorem 
(4.  1)  hold. 

B 

Proof:  See  Appendix. 


V.  Conclusion 

A  class  of  positive  real  multipliers  is 
obtained  to  establish  frequency  domain  condi¬ 
tions  for  stability  of  feedback  systems 
containing  magnetic  hysteresis  non-linearity. 
To  obtain  the  results,  the  model  of  Chua  and 
Stromsmoe  [11]  for  hysteresis  is  employed 
and  the  concepts  of  positivity  and  passivity  are 
utilized. 


VI.  Appendix 

To  9imp\iSy  the  proof  of  Theorem  (4.  1)  and 
Corollary  (4.  1),  the  following  two  lemmas  will 
be  proved  first. 

Lemma  (A.  1):  Let  g€  sector  ((),•), 

>  0.  tod  Men.  Thai  th.  system  of 
Fig.  (A.  1)  Is  pe. sire. 
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Fig.  (A.  1) 


Proof:  For  *11  x€L2e  and  T  >  0 
y>T  =  <x,gM*1x>T.  Let 

M‘*x  =  x  , 

=  <Mx,  gx>j 

=  <  eoS  X.  gx>T  +  <  eorQx.  gx  >T 

♦  <  £(S)x,  gx  >j  (A .  1 ) 

where  Sx  =  ^{x)#  gx  =  g(x),  end  I(S)  is  the 

rational  part  of  M(a).  The  first  two  inner- 
products  of  equation  (A.  1)  are  passive  [4]. 

<£(S)x,  gx  >T  =  <  PTI(S)P_x.  P  gP  x  > 

*  <  PTx,E*(S)P^gPTx> 

=  <  PTx.r*(S)gPTx> 

But  P^.x€L2.  thus,  I*(S)  is  a  mapping  of  L2*-L2. 
Therefore,  by  Theorem  (3.  2) 

<  PTx,  I*(S)g  PTx  >  >0 

which  implies  passivity  of  Fig,  (A.  1).  g 

Lemma  /.A.  2):  Let  f(  •)€(<).  »),  d^x’  >  0, 

and  Then  the  system  of  Fig.  fA,  2)  is 

passive.  g 


But  P^S 


*  PT*  X  SXj. 


TSxT 

=  <  Siij.,  US)  l  xT  > 

*  <  E*(S)  Si,..  fxT> 

Substitution  of  E*(S)  Implies  that:  SZ*(S)  -  I  or.S 
•  »  B  Sl(P;**nK'S| 

♦  Z  (S)  where  I*(S)  §  I  —  i  -° -  .  Then 

i.l  -S"Pi 

*<J1aisZr-£Zr>  +  <  > 


The  first  inner -product  is  passive  [4], 
while  the  second  one  is  positive,  by  Theorem 
(3.  2).  Thus  the  passivity  of  Fig.  (A.  2)  follows. 


Fig.  (A.  3) 

Proof  of  Theorem  (4.  1):  By  inclusion  of  the 
multiplier  M(s),  transform  the  feedback  system 
of  Fig.  (2)  to  the  one  shown  in  Fig.  (A.  3).  For 
the  feed  forward  block: 


Fig.  (A.  2) 


Proof:  For  all  x(t)€L2  and  T  >  0, 

<  x,y  >T  =  <  x,  Mf  ^  x  >T  where  ^  x  =  J x  dt. 

Let  jx  =  x,  then  x€ L2e  since  x€  L2ft. 

*  <Sx,  Mfx 

=  <  Sx.  B  Sfx  >T  +  <  Sx,  6  :  fx  > 

0  1  O  O  1 

+  <  Sx.  E(S)£x  >T  (A.  2) 

The  first  two  Inner -products  of  equation  (A. Z) 
are  passive  [4]. 

<  Sx,  E(S)£x  >T  x  <  x.  E(S)  £  j  x  >T 

*  <  PTx.  PTI(S)f  5  PTx  > 

L«t  ij.  *  g-  P^x,  than  Sx.j.  *  P^.x  €  L 2 

•  <35^.,  PjE(S)fx- > 

*  <PTS  E(S)f  xT  > 


<y1.‘1  >T  =  <yr;>T  +  <yj.?>T 

By  Lemmas  (A.  1)  and  (A.  2),  the  inner  products 
on  the  right-hand  side  of  the  above  equality  are 
passive.  Thus,  the  feed  forward  block  of 
Fig.  (A.  3)  is  passive. 

The  feedback  block  of  the  above  system  is 
strictly  passive  with  finite  gain  if  the  condition? 
set  on  H(s)  and  inequality  (4.  1)  are  satisfied. 

Since  Uj  =  m(t)  *  Uj,  where  m(t)  is  Inverse 
Laplac e Transform  of  M(s),  and  u^,  €L2,  thus 

u1(t)€L2>  Therefore,  by  Theorem  (3.  1).  the 

feedback  system  of  Fig.  (A .^3)  is  finite  gain 
L2-stable,  i.  e.  ,  e^.  e2*  y  ^  y2€  L2  and 

!l«ill.  tl«2(i.  Ily|fi.  Ily2ll  £  p’1fl“1(t>H  +  p2I!u2h 

But 

||S1(t>|i  =  ||  m(t)  *  u j  (t) || 

<  p1II'*1II  +  Pjllujl! 

The  last  inequality  is  obtained  by  substitvxion 
of  m(t>  and  simplification  of  the  norm.  Thus: 

l«tl.  \\*2\\.  lly,B.  Ilf2fl <  pjll«,'l+  pjl'i, 'i.  d2"u2'i 


ib 
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From  Fig.  (A.  3),  y^s  m~*(t)  m*^(t)€Lj 

and  ?2^L2  implie»  that  y2(t),  y2(t)€L2  (*2» 

Appendix  C].  Furthermore,  llyjll^  -  ]|m~*(t)* 

<  ||m"1(t)H  L  ||y2(t)|]L  [12,2App.C]. 

1  2  ^  .1  ^ 

m'1  (t)€L.  implies  that  ||m  (t)]|t  =  C,  a  constant. 

Thu.  1  **» 

l|y2(t)!!L2<  Ci|y2llL2 


=  Cpjllujll+Cp’jllujll  +  Cpjiluj!! 

On  the  other  hand,  y 2=  m*^(t)*y2.  Similar 
upper  bounds  for  *re  imme^i*te,  since 

m'^tKL.. 

1  -1  - 

Similarly,  e^m  (t)*ej.  Thus,  similar 
conclusions  as  y2  for  e^  follow  immediately. 


b)  y 2*  1*2  and  .  e  €L2  imply  that 

y.  e€Lw.  are  continuous,  and  go  to  zero  as  t  ~  • 

[13].  Since  the  model,  i.  e.  ,  equation  (2.  1).  is 
a  continuous  mapping  from  input  to  output  [21]. 
therefore,  e^L#,  and  e^t)-  0  as  t  -  •  imply 

that  the  same  properties  hold  for  y^(t),  i.  e.  , 
yj(t)€Lwf  is  continuous,  and  go  to  zero  as  t 
Similar  conclusions  for  e^  are  immediate. 

■ 


Proof  Outline  of  Corollary  (4.  1);  Apply  a 
positive  feedback  of  gain  j  around  g.  To 
compensate  for  it.  apply  a  positive  feed¬ 
forward  with  gain  £  to  H<s).  Let  g  =  (g~*- 

Then  g€  sector  (0,  ®),  is  monotone  increasing, 
and  g(0)  =  0.  Follow  the  same  procedure  as 
Theorem  (4.  1).  Conclusions  are  immediate. 
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